The Nature of the Condensate in Mass Transport Models 
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We study the phenomenon of real space condensation in the steady state of a class of one 
dimensional mass transport models. We derive the criterion for the occurrence of a condensation 
transition and analyse the precise nature of the shape and the size of the condensate in the condensed 
phase. We find two distinct condensate regimes: one where the condensate is gaussian distributed 
and the particle number fluctuations scale normally as L 1 / 2 where L is the system size, and a second 
regime where the particle number fluctuations become anomalously large and the condensate peak 
is non-gaussian. We interpret these results within the framework of sums of random variables. 
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Condensation transitions are ubiquitous in nature. For 
systems in thermal equilibrium, clustering is well under- 
stood in terms of the competition between entropy and 
energy (typically associated with attractive interactions). 
More exotic are condensations in systems with no in- 
teractions, e.g., Bose- Einstein's free (quantum) particles. 
Less understood are such transitions in non- equilibrium 
systems, in some of which even the concept of energy is 
dubious. For example, condensation is known to occur 
in many mass transport models, defined only by a set 
of rules of evolution, with no clear 'attraction' between 
the masses [1-6]. The relevance of these models lies in 
their applicability to a broad variety of phenomena, e.g., 
traffic flow [7] force propagation through granular media 
[9], granular flow [10] and network dynamics [11]. Cor- 
respondingly, the condensation transition describes jam- 
ming in traffic [3] , bunching of buses [5] , clogging in pipes 
[5], coalescence of shaken steel balls [10] and condensa- 
tion of edges in networks [11]. 

How such transitions arise is especially intriguing for 
one dimensional (d — 1) systems with local dynamical 
rules. A well known example is the Zero-Range Process 
(ZRP) [2,12,13] in which masses hop from site to (the 
next) site according to some transfer rule. In the steady 
state, a finite fraction of the total mass 'condenses' onto 
a single site when p, the global mass density, is increased 
beyond a certain critical value: p c . The system goes from 
a fluid phase, where the mass at each site hovers around 
p, to a condensed phase, where a fluid of density p c co- 
exists with a condensate containing all the 'excess' mass. 

Though condensation in these systems share interest- 
ing analogies [3,6] with the traditional Bose-Einstein con- 
densation, there are important differences. For example, 
here condensation occurs in real space and in all dimen- 
sions. Moreover these systems are non- equilibrium in the 
sense that they are defined by the dynamics, generally 
lack a Hamiltonian and the stationary state is not spec- 
ified by the usual Gibbs-Boltzmann distribution. There 



are two major problems that one faces in the analysis of 
condensation phenomenon in these systems. First, the 
stationary state itself often is very difficult to determine 
and secondly, even if it is known such as in ZRP, the anal- 
ysis of condensation has so far been possible only within 
a grand canonical enemble (GCE) where one is already 
in the thermodynamic (L — > oo) limit. While the GCE 
approach correctly predicts when a condensation transi- 
tion can happen and even the value of the critical density 
Pc it fails to provide much insight into the 'condensed' 
phase (p > p c ) itself. For that one needs to work in a 
canonical ensemble with the system size L finite, which 
has not been possible so far. In this Letter, we show that 
both of these problems can be overcome in a general class 
of mass transport models recently introduced by us [14]. 
This allows us to explore the condensed phase in detail 
revealing rather rich physical behaviors, in particular the 
existence of two different types of condensates. 

Our model is defined as follows: a mass m, resides at 
each site i of a d = 1 periodic lattice of size L. At each 
time step, a portion, m t < m;, chosen from a distribution 
(j)(fh\rn), is chipped off to site i + 1. The dynamics con- 



serves the total mass M = ^2% 



= pL. The model is 



general enough to include many previously studied mod- 
els as special cases [14]. Choosing the chipping kernel 
<p(m\m) appropriately, recovers ZRP, the Asymmetric 
Random Average Process [15] and the chipping model 
of [6]. Moreover the model encompasses both discrete 
and continuous time dynamics and discrete and continu- 
ous mass. In particular, it was shown that the stationary 
state has a simple, factorised form provided the kernel is 
of the form <j>(fh\m) cx u(m)v(m — rfi) [14], where u(z) 
and v[z) are arbitrary non-negative functions. Then the 
joint distribution of mass in the steady state is given by 
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where /(to) = f™ dmu(m)v(m — to) and the 'canonical 
partition function' Z(M, L) is just the normalization 



Z(M,L) = Y[ dmif(mi)S 

1=1 Jo v=i 
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(2) 



Note that (1) is a product of single-site weights f{rrii) 
but the (^-function in (1) implies a fixed total mass thus 
inducing correlations between sites and in general the 
single-site mass probability distribution p(m) ^ f(m). 

The dynamics of the model specifies the functions u(m) 
and v(m — rh), which in turn specify the steady state 
uniquely in terms of weight function f(m). Having de- 
termined the steady state, one next turns to the issue 
of condensation. In particular, we ask: (i) when does a 
condensation transition occur (ii) if condensation occurs, 
what is the precise nature of the condensate? 

The factorization property allows (i) to be addressed 
rather easily within a GCE framework - a la Bose Ein- 
stein. The approach implies taking the L — > oo limit 
and setting the single-site mass distribution function 
p(m) — /(m)e~ Mm where p is the chemical potential and 
is chosen to fix the density p = J dnip{m)m. Thus, con- 
densation must occur for p > p c ~ J dmf(m)m which is 
the maximum allowed value of p within the GCE. Based 
on previous works on the ZRP related case [4,5], it is 
easy to show that a condensation transition occurs if the 
single site weights decay for large to as 



/(to) ~ Am 7 with 7 > 2 



(3) 



A simple example of a chipping kernel which gives such 
weights is furnished by u(fh) = exp(— am) and v(m — 
to) = (1 + to — m)~ 7 which yield f(m) ~ to ' /a for large 
to. In the following, we stay with the choice of /(to) in 
(3) and set, without loss of generality, J °° f(m)dm = 1. 

The GCE analysis correctly predicts the criterion for 
condensation and even the critical density p c , but pro- 
vides little insight into the condensed phase itself where 
p > p c . In this work, we are able to explore the 
condensed phase by staying within the canonical en- 
semble and analyzing the mass distribution p(m) = 

J dm2----dm L P(m, • • • , m L )S (j2j=2 m j + 77 
finite system of size L. Using (2), we have 
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p(m) — /(to) 



Z(M — to, L — 1) 
Z(M, L) 



(4) 



The rest of the letter is devoted to the analysis of p(m) 
in (4) with /(to) given by (3). We have two parameters 7 
and p. Our goal is to show how the condensation is man- 
ifested by different behaviors of p(m) in different regions 
of the {p — 7) plane giving rise to a rich phase diagram 
in Fig. 2. 

First, consider the Laplace transform of (2): 



/•OO 

/ Z{M,L)e- sM dM =[g{s)] L , 
Jo 



(5) 



where g(s) = J °° f(m)e~ sm dm. The main challenge is 
to invert (5) for a given /(to) and exploit its behavior 
to analyse p(m). Before proceeding to the general case, 
let us present a case in which both Z(M, L) and p(m) 
can be obtained in closed form. Wc choose /(to) = 
2e- 1 A™TO- 5 /70F, for which ff ( s ) = (1 + 2s 1 / 2 )c- 2s1 ' 2 
and our results below show that p c = —g'(0) = 2. Now, 
(5) can be inverted to provide the closed form [16]: 



(6) 



Z(M, L) - B m ,l H L (r) - VM w (r 



where B M , L = LAf-( i+3 )/ 2 e~ L l M /y^r, r = (2L + 
M)/2y/~M and Hk{r) is the Hermite polynomial of degree 
k. Substituting Z(M, L) in (4), we obtain p(m) explicitly 
and plot the case of L = 100 in Fig. 1. The transition 
from the subcritical (p = 1), through the critical (p = 2), 
to the supercritical (p — 6) cases, is clearly seen: the con- 
densate showing up as an additional, asymmetric bump 
for p = 6. The explicit solution in this toy example pro- 
vides us with useful insights into what to expect in the 
general case. 







L=100 




p=2 






p=l \ 








P <> i 



In (m) 

FIG. 1. The distribution p(m) vs. m for the exactly solv- 
able case, plotted using Mathematica for L — 100 and p = 1 
(subcrtical), p — p c — 2 (critical) and p = 6 (supercritical). 
The condensate shows up as an additional bump near the tail 
of p(m) in the supercritical case. 

Before proceeding, let us summarise here our main re- 
sults: we refer to Fig. 1 for typical forms of the mass dis- 
tribution p(m) and Fig. 2 for a schematic phase diagram. 
In the subcritical regime the system is in a fluid phase 
where the mass distribution decays exponentially with 
decay length increasing with density. At p c the distribu- 
tion decays as a power law p{m) ~ m~ 7 and at p > p c 
the distribution develops an extra piece, representing the 
condensate, centred around M — Lp c . By our analysis 
within the canonical ensemble we show that this piece 
will have a 'normal', gaussian form when 7 > 3. When 
7 < 3, however the condensate will have an anomalous, 
asymmetric form as seen in Fig. 1 for p = 6. In the 
following we supply explicit expressions for these forms. 
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FIG. 2. Schematic phase diagram in the p-7 plane, 
full line represents the critical density p c ("i) 
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We formally invert (5) using the Bromwich formula, 



Z(M, L) 



s +ioo 



ds 



,exp[L(\ng{s)+ps)} (7) 



where the contour parallels the imaginary axis with its 
real part, s , to the right of all singularities of the inte- 
grand. Since /(to < 0) = 0, the integrand is analytic in 
the right half plane. Therefore, so can assume any non- 
negative value. Meanwhile, for / given by (3), s = is 
a branch point singularity. As we shall see, in the sub- 
critical case there exists a saddle point at positive s and 
so can be chosen to be this saddle point, whereas in the 
critical and supercritical cases the leading contribution is 
obtained by wrapping the contour around s = 0. 

First we evaluate (7) in the limit L — > 00 by the saddle 
point method, assuming it exists. Let h(s) = ps + lng(s). 
Then the saddle point equation, h'(s ) = 0, is 



-g'(s )/g(s ) = pOo) 



leading us to, e.g., 



Z(M,L) ~ 



exp(£/i(s )) 

^/2TTLh"( SQ ) ' 



(8) 



(9) 



If p < p c = p(0), then (8) has a solution for s > 0, 
the saddle point approximation is valid, no condensation 
occurs. Substituting (9) in (4) we get, for p < p c and 
to <C (/9 C — p)L, p(m) ~ f{m)c~ s ° m , recovering the GCE 
upon identifying the chemical potential p — s n . 

We now focus on the behavior as we approach critical- 
ity from the subcritical regime and consider (8) for small 
and positive p c — p, i.e., small sq. Thus, we just need the 
small s behavior of h(s). For /(to) in (3) with a non- 
integer 7, one can expand, quite generally, the Laplace 
transform g(s) of (3) for small s, as 



n-l 



fe=0 



(10) 



Here n = intfy], pt is the k th moment of /(to) (which 
exists for k < n). The second term of (10) is the leading 



singular part and it can be shown that b = AT (I — 7). 
Note that /zo = 1 for normalized /, p\ — —g'(0) = p c 
and A = \J P2 — Mi, which is the width of the distribu- 
tion /, is finite if 7 > 3. The role of 7 = 3 is now clear. 
The next-to- leading term is s 2 in one case and s 7_1 in 
the other, so that the saddle-point solution in (8) is given 



(Pc-P) 



6(7-1) 



1/(7-2) 



by, s ~ (p c - p)/A 2 for 7 > 3 and s 
otherwise. Inserting this behavior in the expression (9) 
gives for 7 > 3, Z(M, L) - exp(-L(p c - pf /2A 2 ), point- 
ing to a system with Gaussian distributions and normal 
fluctuations. In contrast, for 2 < 7 < 3, wc will show 
that anomalous fluctuations and non-Gaussians appear. 

For the supercritical regime (p > p c ), there is no so- 
lution to (8) on the positive real axis and more care is 
needed to find the asymptotic form of Z(M, L). Our ap- 
proach is to use the fact that the integral (7) will be dom- 
inated by s ~ 0. Thus we can use the small s expansion 
(10) and develop a scaling analysis by identifying the dif- 
ferent scaling regimes and calculating the corresponding 
scaling forms for Z(M, L) in the large L limit. 

Using (10), one can rewrite (4) as 



p(m) ~ /(to) 



W({m-M ex )/L) 
W (p c - p) ■ 



(11) 



where 



W(y) 



J -too 2tti 



exp 



L{-y S +^-s 2 + ... + bs^- 1 ) 



and M ex = (p — p c )L is the excess mass. All crucial 
information about the condensate 'bump' is encoded in 
the asymptotic behavior of W. Deferring the details to a 
later publication [16], we outline our main results below. 
Again, we consider the two cases (7 > or < 3) separately. 

Case-I (7 > 3): We find that, for large L, and in the 
Oil}! 2 ) neighborhood of M ex , the condensate appears in 
p(m) as a pure Gaussian and can be cast in scaling form: 



Pcond(TO) 



1 



-z 2 /2. 



z = 



AL 1 / 2 



(12) 



Note that its integral over to is 1/L, indicating that the 
condensation occurs at a single site. 

Case-II (2 < 7 < 3): In the neighborhood of the excess 
mass M ex , we find that p(m) has the scaling form: 



M, 



£1/(7-1) 



(13) 



where VJz) = f°° ^ 

' v ' J —tOO 27TI 



-zs+bs 1 



Though we have no 



closed form for V 1 {z), we obtain its asymptotics: 



Vy(z) ~ as z 



-00 



ClZ (3-7)/2(7-2) g-o: 



(T-l)/(T-2) 



(14) 

as z — > 00 (15) 



where c\, ci are constants dependent on 7. Note that 
this condensate 'bump' is far from being gaussian: it has 



3 



a highly asymmetric shape, evidenced by (14, 15). The 
peak occurs at m = M ex and scales as ~ £-7/(7-1) 
Meanwhile its width is 7 dependent: L 1 ^ 7-1 ). The area 
under the bump is 1/L, again implying that the conden- 
sate occurs at only one site. 

Finally, in both cases I, II the supercritical partition 
function is given by 

Z(M,L)~AL/Ml (16) 

Similar results for the critical case p = p c and 7 = 2, 3 
(where one obtains logarithmic corrections) will be pub- 
lished elsewhere [16]. 

The implication of these results are clear: In the con- 
densed phase, the condensate acts as a reservoir for the 
critical fluid. Thus, the width of the condensate bump re- 
flects the mass fluctuation in the critical fluid. For 7 > 3 
we showed that condensate is gaussian distributed with 
width AL 1 / 2 , so that the masses in the fluid fluctuate 
normally. For 7 < 3, however, the width of the con- 
densate and the mass fluctuations in the fluid are both 
anomalously large, namely, 0(L 1 ^ 1 ~ 1 ' > ). Further impli- 
cations concern the dynamics within the steady state. In 
systems with symmetry breaking, the 'flip time' r [17] is 
of interest. Here the translational symmetry is broken 
by the selection of a site to hold the condensate and r 
corresponds to the typical time a condensate exists be- 
fore dissolving and reforming on another site. A rough 
estimate for r is Pco n d( m )i with \iti — M ex \ <~ O(L). For 
7 > 3, this implies flip times growing exponentially with 
the system size, whereas for 7 < 3, they would diverge 
more slowly, as some power of L. 

Our results may be naturally interpreted within the 
framework of sums of random variables. The parti- 
tion function (2) is proportional to the probability that 
the sum of L independent random variables m„ each 
distributed according to f(m) is equal to M. Given 
/(to) ~ AmT 1 for large to, the rrii are thus Levy flights 
and Z(M, L) is just proportional to the probability dis- 
tribution of the position of a Levy walker (only taking 
positive steps) after L steps. Then (16) can be inter- 
preted in terms of the extreme statistics of a Levy walk. 
The mean of a sum of L random variables is just p,\L. 
Thus if M < Lpi = M c one expects the sum to con- 
tain random variables of typical size O(l), whereas for 
M > M c one expects the sum to be dominated by a 
rare event i.e. L — 1 of the variables would be of order 
O(pi) except for one which would be large and equal 
to M — M c . Given that the distribution of the random 
variables is <~ AmT 1 ', the probability that this large vari- 
able takes the value M - M c is A(M - M c )~ 7 . This 
large contribution could be any of the L possible ones, 
thus the total probability is AL(M — M c )~ 7 , recovering 
(16). Moreover, the gaussian L 1 / 2 fluctuation for 7 > 3 
and non-gaussian L 1 /^-!) fluctuation for 2 < 7 < 3 cor- 
respond respectively to the normal and the anomalous 
diffusion of a Levy process. 



To summarise we have considered a very broad class 
of mass transport models and derived the condition for 
condensation. We have presented an analysis within the 
canonical ensemble that elucidates the nature and struc- 
ture of the condensate. In particular we have identified 
two distinct condensate regimes where the condensate is 
normal and anomalous and derived the scaling distribu- 
tion for the two types of condensate. Our results rely on 
the factorisation property of the steady state (1), but we 
believe the phase scenario of Fig. 2 may apply in models 
without factorised steady states; it would be of interest 
to verify this. The underlying dynamics of the model 
we have studied are one-dimensional however the condi- 
tion for a factorised state may be generalized to higher 
dimensions [16] and work is in progress to investigate 
condensation in such systems. Our results confirm that 
condensation may occur in a wider class of continuous 
mass models as well as discrete mass models such as the 
ZRP, as suggested in [18]. It would also be of interest to 
use this class of models to generalize the phase separation 
criterion of [8] which is based on the ZRP condensation. 

This research is supported in part by the US NSF 
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